Hydrodynamic description of Hard-core Bosons on a Galileo ramp 
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We study the quantum evolution of a cloud of hard-core bosons loaded on a one-dimensional optical lattice 
after its sudden release from a harmonic trap. Just after the trap has been removed, a linear ramp potential 
is applied, mimicking the so called Galileo ramp experiment. The non-equilibrium expansion of the bosonic 
cloud is elucidated through a hydrodynamic al description which is compared to the exact numerical evolution 
obtained by exact diagonalization on finite lattice sizes. The system is found to exhibit a rich behavior showing 
in particular Bloch oscillations of a self-trapped condensate and an ejected particle density leading to two di- 
verging entangled condensates. Depending on the initial density of the gas different regimes of Josephson-like 
oscillations are observed. At low densities, the trapped part of the cloud is in a superfluid phase that oscillates 
in time as a whole. At higher densities, the trapped condensate is in a mixed superfluid-Mott phase that show a 
breathing regime for steep enough potential ramps. 
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I. INTRODUCTION 

It is of primary importance to understand the dynamical be- 
havior of strongly correlated many-body quantum systems, 
especially since their behavior in many cases could show 
counter-intuitive effects with respect to the classical naive ex- 
pectation. One of the early discovered such effect is the so- 
called Bloch oscillation regime occurring when a small con- 
stant force F is applied to a quasi-particle living on a lattice 
dl. In such a case, the particle momentum is drifted in time 
according to q(t) = ^(0) + Ft modulo the Brillouin zone, 
leading to a periodic motion with period tb = 27r/\F\ (h and 
the lattice spacing are set to one). Such oscillatory behav- 
ior has been observed in many different physical contexts like 
semiconductors, thermal gases, photonics, cold atoms physics 
and with phonons |2J. At stronger forces, in multi-band sys- 
tems, such oscillations can be suppressed by the Landau- 
Zener interband tunneling O. The Bloch oscillations can also 
be suppressed due to many-body effects with damping related 
to the integrability of the system (HISl. 

On the other hand, experimental advances, notably on ul- 
tracold atomic gases have lead to a revival of the field of low- 
dimensional many -body systems |6|, especially within non- 
equilibrium aspects in in. Thanks to the very low dissipa- 
tion rate and very long time phase coherence of such sys- 
tems, it has become possible to realize experimentally almost 
the unitary dynamical evolution of ideal models like the one- 
dimensional Bose-Hubbard model, even in its hard-core limit 
ll6ll . In those systems, the release of the gas from a trap, 
which is a standard procedure for cold atoms experiments, 
can lead to interesting metastable states 1 9 1 or even to dynam- 
ical fermionization effects when the frequency of a parabolic 
trap is modified fTO]. Recent studies have considered the 
dynamical behavior of hard-core bosons in tilted optical lat- 
tices HTHIJI, focusing in particular on the Bloch oscillations, 
Landau-Zener tunneling and entanglement properties. 
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In this study we pursue on the direction set by 1 12] by con- 
sidering now a hard-core boson cloud initially localized on 
part of an optical lattice thanks to a harmonic trap. The trap is 
then suddenly removed and a linear ramp potential (a constant 
force F) is added to the region initially occupied by the cloud. 
The unitary expansion of the cloud leads to the escape of part 
of the condensate while the remaining part is self-trapped into 
the initial region and performs Bloch oscillations. Thanks to 
the initial correlations within the cloud the escaping particles 
are strongly entangled with the localized oscillating conden- 
sate. Depending on the initial density, two different situations 
arise. Namely, at low density the trapped condensate is in 
a superfluid phase that oscillates in time as a whole while 
at higher density it is in a mixed superfluid-Mott phase that 
show a richer behavior with, in particular, a breathing regime 
at large enough forces. The out-of-equilibrium dynamical be- 
havior of the system is analyzed numerically by means of ex- 
act diagonalization following the methods developed in |[T3l 
and through an analytical hydrodynamical approximation al- 
ready used in |[T2l . 

The paper is organized as follows: in section II, we present 
the model, its mapping to a Fermi system and the hydrody- 
namical description used in this paper. In section III, we study 
the dynamics after the sudden release of the bosonic cloud 
focusing first on the emitted wave-packet, then on the self- 
trapped condensate and finally discussing the entanglement 
between the emitted particles and the self-trapped ones. In 
section IV, a brief summary is given. 



II. THE MODEL AND ITS HYDRODYNAMIC 
DESCRIPTION 

A. The model 

We consider in this study the one-dimensional Bose- 
Hubbard model describing a set of bosons leaving on a lattice 
with a repulsive on-site interaction U and submitted to an ex- 
ternal potential V{t) which may vary in time. The dynamics 
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is generated by the Hamiltonian 



U 



H = -J^[a]+ia,+h.c.] + -^n,(n,-l)+^F,(t)n 



with aj and the usual destruction and creation bosonic op- 
erators, rij = a^-aj the bosonic density at site j. J is the hop- 
ping magnitude and will be set to J = 1 /2 in the rest of this 
work (this will give a band width A = 1). In particular, we 
focus our attention on the hard-core boson limit of this model 
by tacking the limit U ^ oo. In this limit, the dynamics can 
be described by the new hamiltonian 



n 



(1) 



where the new bosonic creation and annihilation operators 
and hj satisfy the on-site anticommutation rule {hj^h^-} = 1, 

{hj^hj} = {6^, 6^} = 0, avoiding a double occupancy of the 
same site, while commuting at distinct sites. The standard 
procedure to diagonalize this Hamiltonian is to fermionise 
it through a Jordan- Wigner mapping and then performing a 
canonical Bogoliubov transformation to new diagonal Fermi 
operators. Indeed, introducing the lattice Fermi creation op- 
erators = ni<j(l ~ 2ni)6^ and their adjoint annihilation 

operators Cj = (c]) + , rejecting the boundaries at infinity, the 
Hamiltonian is expressed as a tight-binding Fermi system 



n 



h\hj the occupation operator at site j. 
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Thanks to the quadratic form of this Hamiltonian, it is readily 
diagonalized through a standard canonical transformation |T3l 
leading to the diagonal expression T-L = "^q^qvlvq where the 
Cq are the excitation energies associated to the free Fermi par- 
ticles created by t^J and destroyed by r]q and the problem is 
in principle solved. However, for a quite general inhomoge- 
neous potential Vj{t) one has to compute numerically on a 
finite lattice the single-particle spectrum e^, while the relevant 
observables are expressed as determinants of the single par- 
ticle Green functions (cjcj) which have also to be computed 
numerically in the general case. Nevertheless, in the limit of 
a very large system with a sufficiently smooth potential func- 
tion V{x,t) one can well capture the features of the boson 
density dynamics by a hydrodynamic limit using a continuous 
description of the model Iil2il . 



B. Continuum limit and local equilibrium hypothesis 

Consider the one-dimensional Hard-core bosons system on 
an infinite one-dimensional lattice with lattice spacing a <C 1. 
Let the potential V{x) he 3. smooth real function (at least let 
us say V G C^) on the lattice. We split the real line into 
regular intervals [x^x -\- Ax], with Ax = aM, containing a 



large number M of lattice sites, while keeping the width Ax 
small enough in the sense that Vja G [x, x + Ax], V{ja) 
V{x), that is the potential keeps almost a constant value on 
each interval. The Hamiltonian can be recast in the continuum 
limit a ^ 0, Ax while keeping Ax /a = M 1 in the 
following form 



/oo 
dx n{x) 
-oo 

where the hamiltonian density 1-L{x) is given by 
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n{x) 
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dy 



^^■^(x + y)'^{x + - a) + h.c. 



^^■'■(x + y)V{x)'^{x + y) 



(4) 



in terms of the continuous creation and annihilation Fermi 
field operators ^^(^) and ^(^). To achieve ^ with ^ we 
have supposed that the interaction contribution between dif- 
ferent intervals, which is just a local boundary term, is very 
small compared to the contribution ^ within the interval. 
This is true if the potential variations are sufficiently small. 
The local Hamiltonian density 1-L{x) is simply a continuous 
version of ^ with a constant potential V{x), and it can be 
consequently easily diagonalized through the canonical map- 
ping 



^(x 



y)= dq (t)q{x^y)r]{q,x) , 

JO 

^\x^y)= / dq(l)l{x^y)r]\q,x) , 
Jo 



(5) 
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where the field operators r]{q^x) and r]^{q^x) annihilates 
and respectively creates in region [x, x + Ax] a particle 
with momentum q and satisfy the anti-commutation rules 
{r]{Pj x)^r]\p' ^ x')} = S{p — p')S{x — x'). The exact form of 
the Bogoliubov functions (l)q{u) entering into the definition of 
the new fields essentially depends on the boundary conditions 
imposed on the interval [x, x + Ax]. Using these new fields, 
the total Hamiltonian takes the diagonal form 



dx I dq \y{x) 



q]r]^q,x)r]{q,x) . (7) 



Therefore, with respect to this Hamiltonian, the -particles 
ground state of the bosonic system is given by a local equilib- 
rium state. That is, all the quasi-particles associated to each 
phase-space points (g,x), with energies e(x,g) = V{x) — 
cos q bellow the Fermi level cpiN), Sive added to the vacuum 
state: 



i^o)= n vHQ.m- 

q,x 



(8) 



The Fermi energy is given by imposing the constraint 
/ dx p{x) = J dx J dq r]^ {q,x)r]{q^x) = N on the total 
number of particles. This readily implies for the ground-state 
density profile 

V(x) - > 1 
;^arccos(l/(x) -ej^) |V(x)-eF|<l (9) 

1 V(x) - eF < -1 



p(x) 
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In Figure ([T]) we compare the exact numerical diagonaliza- 
tion with the local equilibrium prediction on a system with 
L = 400 sites at half filling and, as an illustration, a poten- 
tial V(x) = 2x/L + 1.5sin(8x/L) + sin(16x/L). Notice the 
good matching between the continuum prediction ^ and the 
numerical data. In particular, we have also graphically repro- 
duced the occupied energy levels (yellow region) at half filling 
factor. In practice, since cosq G [—1,1], the energies e(g, x) 
fall in the gray strip which follows the shape of the potential. 
The local equilibrium approximation breaks down whenever 
the local number of particle is small and the potential varies 
sharply. 




FIG. 1. (Color online) Ground state density profile for a system 
with a potential V{x) = 2x/L + 1.5sin(8x/L) + sin(16x/L) and 
N = 200 (Exact numerical diagonalization results are represented 
by the symbols while the local equilibrium distribution is given by 
the dashed line). The local energy band is represented together with 
the local one-particle filled states up to the Fermi level e n • 

In a grand canonical situation, when the number of particles 
is not fixed, the Fermi level ep in the ground state as to be set 
to since the ground state is build up by adding to the vacuum 
all excitations with negative energies : 



i^o)= n vHq.^. 

{iq,x)\eiq,x)<0} 



(10) 



and consequently the density profile is given by ^ with 
ep = and the total number of particles in that state is just 
Jdx p{x). 



with A < and (Iq > —1 such that there is at least few par- 
ticles loaded on the lattice. This parametrization insures that 
the bosonic condensate is centered at A/2 with a spacial ex- 
tension of width \A\ such that outside the region x G [A, 0] 
the density (in the hydrodynamic limit) is vanishing (see fig- 
ure [2]). For a given width | A| of the condensate, we have two 
qualitative distinct situations that are controlled by the value 
of the chemical potential /io- Indeed, for jiQ < 1, the poten- 
tial is never smaller than —1 and the condensate is fully in its 
superfluid phase with a density profile 
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(12) 

where 6{x) is the Heaviside function. We will refer to this 
state as the pure superfluide phase (SF-phase). On the other 
hand, for /io > 1, there is in the middle of the condensate a 
Mott insulating phase, with p = 1, that extends spacially over 
the region x G [f - Amou, f + ^Mott] with 
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(13) 



This Mott phase is surrounded by two superfluid phases, of 
spatial extensions 
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with a profile given by ([9| with the harmonic potential ( pT] ), as 
illustrated on figure [2] As the chemical potential /io is getting 
larger and larger, the Mott phase is growing at the expense of 
the superfluid phases, shrinking them to the boundaries close 
io X = A and x = 0. We call this state the mixed superfluid- 
Mott phase (Mixed SF-Mott phase). 




FIG. 2. (Color online) Left: Density profile of the initial SF state 
with \A\ — 200 and /io = 0. Right: Density profile of the initial 
SF-M state with \A\ = 200 and /io = 3. The symbols are obtained 
by exact diagonalization while the full lines give the hydrodynamical 
profile. 



C. Initial state 



In the following we will fix the initial state of the bosonic 
condensate as the grand canonical ground state associated to 
a harmonic trap potential. The harmonic potential is parame- 
terized as 



V{x,t< 0) = a 



/io , a 



A^ 



(l + Mo) , (11) 



III. DYNAMICS AFTER THE SUDDEN QUENCH 

A. Sudden release of the trap and loading of the linear ramp 

At time t = 0+, starting from the previous initial state, we 
suddenly release the gas from the parabolic trap and load an 
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external constant force F on the negative side of the real axis. 
The force is described by the linear ramp potential 

V{x,t > 0) = -Fx 0{-x) 

and the system dynamics is governed by a new Hamiltonian 
^ with potential V{x^ t > 0) = —Fx 0{—x). As we will see 
in the following, the main features of the dynamics are well 
described in the hydrodynamic limit. As seen previously, the 
initial state has approximately a coarse-grained phase-space 
density 

w{x, q) = -0{-x)0{x + A)0{q)0{qF{x) - q) 

TT 

with a local Fermi wave- vector qpix) = 7rp(x) where p{x) 
is given by ([9]). The local initial density is then simply 
given by adding all the local quasiparticles (associated to each 
phase-space point) with different wave- vectors q: p{x) = 
J dq w{x,q). Just after the sudden quench, the energy of a 
particle initially at position x is shifted by the additional po- 
tential energy V{x^t > 0). This brutal change of potential 
leads, as a consequence, to a unitary out-of-equilibrium evo- 
lution of the condensate from its initial state. In order to un- 
derstand that evolution, taking into account that the dynamics 
is unitary, one has to realise that the initial quasiparticles that 
build up the local density are emitted to the right and to the 
left on trajectories of constant energies. This approach will 
allow us to reconstruct the evolving density profile p{x; t) and 
the associated current density profile t). 

B. Wave-packet emission 

1. Density profile 

The typical evolution of the bosonic density, just after the 
sudden quench, is shown on figure[3]for several forces and for 
the two distinct typical initial states, that is the pure SF-phase 
and the mixed SF-Mott state. As we see clearly on figure |3] 
only part of the bosons are ejected to the right with a corre- 
sponding density that is spreading in time (notice that since 
the numerical results are obtained on a finite lattice with open 
boundary conditions the ejected particles are reflected when 
they reach the right boundary wall). The remaining particles 
are self-trapped into the initial region, performing Bloch os- 
cillations (we will come back on this point in the next sec- 
tion). The reason for the escape of the particles to the right 
is best understood from the hydrodynamical description. In- 
deed, the population living initially on the tilted band with 
energies within the interval [—1,1] are connected at t = 0+ 
to the propagating states of the energy band at x > (see 
figure |4]a)). Those particles will consequently escape from 
there initial position and propagate to the right toward +oo. 
Right movers will escape directly while left movers will be 
first reflected on the left tilted band edge and then propagate 
toward the right direction. From energy conservation, within 
the tilted band, as for Bloch oscillations, the momentum of a 
right(left) mover with energy e evolves as q^ (t) = ±q -\- Ft, 
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FIG. 3. (Color online) Snapshots of the evolution of the density pro- 
file obtained from exact numerical calculation for the two different 
initial states (SF-Mott and SF phases) at two different forces F. Up- 
left: SF-Mott phase with F = 0.1. Up-right: SF-Mott phase with 
F = 0.04. Down-left: SF phase with F = 0.1. Down-rigth: SF 
phase with F = 0.05. 

describing the conversion of the potential energy V{x) into ki- 
netic energy — cos q(t). As the escaping left and right movers 
have reached the origin x = 0, all the potential energy has 
been converted into kinetic energy. This happens at times 
= {tq-\- arccos(— e)) for the right(+ sign) and the 
left(— sign) movers respectively. In the regi on x > 0, the 
particles evolve at constant velocity v{e) = \/\ — depend- 
ing only on there (initial) energy e. Consequently the initial 
right and left movers are traveling along the trajectories 

x'^{e,q,t) = \/l - e2 {t - t^) (15) 
which can be written in terms of the initial position x as 

x^{x^q^t)= y^l — {Fx + cos^)^ 

X ~ ^ ^ arccos(Fx + cosq))^ . (16) 

The density profile of the escaping particles is then given by 
the sum of the left and right movers contributions 

Pesc(^; t) = t) + p^sci^; t) (17) 

with 

Ptsci^'.t)= f ^1 ^S{x-x^{e,q,t)) (18) 

J-l ^ JQie) 27r 

where the domain of integration Q{e) over the momenta is 
shown on figure [4] a). In figure [4]b) we show the exact numer- 
ical results obtained by exact diagonalization of the escaping 
density compared to the hydrodynamic prediction ([18]). We 
see a very good agreement up to small interference effects 
that are obviously neglected in the continuum limit. It is inter- 
esting to notice that the rightmost front of the escaping wave- 
pacquet shows, beyond the hydrodynamic envelop, a staircase 
structure as already noticed in a different context |14, 151. 
Each of the stairs corresponds to a particle, the integrated den- 
sity over each such stair being equal to one, ejecte d to the right 
and moving ballistically with velocity v{e) = \/\ — e^. Be- 
hind the front the exact structure is a bit more complicated and 
it is difficult to discriminate between different particles. 
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FIG. 4. (Color online) (a) Schematic picture of an initial SF-M state 
just after the sudden quench of the linear ramp. The dashed blue 
region represents the density that will be ejected into the propagative 
band, (b) Density profile of the ejected particles at times, from left 
to right, t = 74, t = 110 and t = 156 . The symbols represent the 
exact diagonalization results obtained on a chain of 400 sites with 
a force F = 0.04 and with an initial SF-M state compared to the 
hydrodynamic predictions (full lines). 



C. Dynamics of the trapped particles 

1. Density profile 

As seen from the snapshot [3] due to energy conservation, 
there is a self-trapped density in the initial region which ex- 
hibits a quite rich behavior with Bloch oscillations related to 
the fact that the particles are traveling on a lattice and experi- 
ence a constant force F 1 1 1. This leads to a momentum drift in 
time according to q{t) = ^(0) + Ft (h and the lattice spacing 
are set to one) modulo the Brillouin zone and then leading to 
a periodic motion with period tb = 2it/F. It is also remark- 
able to notice that while in the case of a pure SF initial state 
the condensate oscillates in time as a whole (both sides of the 
condensate are moving in phase) while in the case of a mixed 
SF-Mott initial state one has to distinguish between a low 
force regime and a high force one. Indeed, at low forces (see 
F = 0.04 in figure [3]) the two superfluid phases surrounding 
the Mott phase are oscillating in phase (both moving on the 
same direction at a given time) while at high forces (F = 0.1 
in figure [3]) the two superfluid parts are moving in opposite 
directions at a given time, leading to a breathing condensate. 
This behavior is easily understood thanks to the hydrodynam- 
ical approach where the density profile of the trapped particles 
is given by the sum of left and right movers densities: 



The left and right movers densities are given by 



(19) 



where 



P^(^'^) = ^ ^^[f2iq)Jliq)]i9^i^^Q^t)) (20) 



(x, q^t) = X — — [cos q — cos{q ± Ft)] (21) 
F 




(22) 
(23) 



with q = arccos(— 1 — Fxm) where is the right-most 
locus of the initial trapped condensate (see figure |4]a)). 

On figure [5] we have plotted over a half-period the oscillat- 
ing density profile of the self-trapped particles for the two ini- 
tial states (SF-Mott and SF phases) as extracted from the exact 
numerical diagonalization results shown in figure |3] and com- 
pared to the hydrodynamical prediction ([19]), together with a 
schematic representation of the initial state. First of all we see 
a very good agreement of the hydrodynamical prediction with 
the exact numerical results. In these figures we clearly see the 
two different oscillation regimes in the SF-Mott case (see fig- 
ures |5]a) and[5]c)). The reason for that is easily understood 
from the schematic pictures of the initial states. Indeed at low 
enough forces, as seen on the schematic picture of the initial 
state, the energy of the local density at the separation point 
xsF-M between the SF phase and the Mott phase is higher 
than the energy at the left-most initial locus x = A and conse- 
quently the density will explore regions on the left of A up to 
the locus XL = xsF-M — '^l F. This phenomenon disappears 
when the energy at the locus xsf-m is getting smaller than 
the energy at the point x = A. This precisely appears above a 
threshold force F* given by -F*A - 1 = -F*xsf-m + 1- 
In the case represented in figure [5]b) the threshold force is 
given by F* = 0.0689. where we see almost no evolution of 
the profile at the left side. Beyond that value, as seen on figure 
|5]c) for F = 0.1, the left and right sides of the density profile 
are moving in phase. With a SF initial state, the self trapped 
density is always globally oscillating as a whole (see figure 

Hid)). 



2. Current density profile 

This remarkable periodic motion of the trapped density is 
naturally associated to a flow of particles giving rise to a pe- 
riodic current density j{x^t) = {J{x^t)) which is defined 
through the continuity equation i[H, n{x)] = — VJ(x) where 
n{x) is the occupation operator at site x. In figure|6]we show a 
snapshot of the current density obtained from exact diagonal- 
ization for a SF state at F = 0.05 (top) and a SF-Mott initial 
state at low forces (F = 0.04) and high forces (F = 0.1). 
In both cases, in the trapped region we see clearly the Bloch 
oscillations with a strip structure for the pure SF phase indi- 
cating the collective motion of the superfluid condensate while 
in the mixed SF-Mott situation we observe the oscillations (in 
phase at F = 0.04 and in opposite phase at F = 0.1) of the 
two SF phases surrounding a stationary Mott plateau. 

In the hydrodynamical limit, the current density of the self- 
trapped condensate is simply given by summing over all the 
quasiparticles current contributions p{q)v{q) with the veloci- 
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FIG. 5. (Color online) (a) Top: Self-trapped density profile for an 
initial SF-Mott state at different times for F = 0.04. The hydro- 
dynamical prediction is given by the dashed line, (b) Middle-up: 
Same as (a) for F* = 0.0689.. (c) Middle-down: Same as (a) for 
F — 0.1. (d) Bottom: Self-trapped density profile for an initial SF 
state at different times for F = 0.05. The hydrodynamical predic- 
tions are represented in each case by the dashed lines. 



ties v{q) — ± sin(g(0) ± Ft) leading to 

dq 



j{x 



^t)= ^sin((7 + Ft)n[j,(^)j^(g)](^+(x,^,t)) 
- ^ ^sin(g-Ft)n[^,(,)j,(,)](^-(x,g,t)).(24) 

We have plotted in figure [7] the current density ( [24] ) compared 
to the exact numerical one at different positions as a function 




FIG. 6. (Color online) Snapshot of the time evolution of the current 
density profile for a SF initial state at F = 0.05 (top) and a SF- 
Mott initial state at low forces F — 0.04 (middle) and high forces 
F = 0.1 (bottom). 



of time for the two cases, SF and SF-M. We observe a per- 



fect matching of the hydrodynamical prediction ( [24] ) with the 
exact numerical values. We also observe a very nice sinus 
behavior of the current in the middle of the superfluid con- 
densate reminiscent of a Josephson-type oscillations which is 
explicitly given from (24) as j(A/2, t) = l/n sm{Ft). 



D. 



Entanglement between the trapped and the escaping 
particles 



Just after the sudden unloading of the parabolic trap and the 
quench of the linear ramp, the initial condensate is split into 
two disjoint parts: the escaping particles and the remaining 
self-trapped ones. Due to the initial correlations in the start- 



ing state (10), these two well separated condensates are en- 
tangled. This entanglement between propagative modes (into 
the right propagative energy band) and bound states (into the 
self-trapping region) can be quantified through the bi-partite 
von Neuman entropy S{x^ t) = — Tr{p(x, t) In t)} where 
t) = TYy^x{\'^{t)){'^{t)\} is the reduced density matrix 
associated to the left of position x on the lattice, deduced from 
the time-evolved pure state by tracing out the degrees 

of freedom on the right of x (see for example the review 1 16 |). 

As seen on figure [8] we have qualitatively two different sit- 
uations depending on whether the initial state is a pure SF 
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FIG. 7. (Color online) Top: Current density at x = —100 for F = 
0.1 (shortest period) and F = 0.05 in the SF case. Bottom: Current 
density at x = -185 for F = 0.1 > F* (shortest period) and F = 
0.04 < F* in the SF-Mott case. Full lines are the corresponding 
hydrodynamical predictions. 




FIG. 8. (Color online) Top: Entanglement entropy as a function of 
space and time for F = 0.3 with a SF initial state generated with 
A = 50 and /xo = 0. Bottom: Entanglement entropy as a function of 
space and time for F = 0.5 with a SF-M initial state generated with 
A = 50 and /xo 2. 



state (see the top snapshot of figure [8]) or a mixed SF-Mott 
one (see the other snapshot of figure jSJ). First of all, in both 
cases we see clearly the entanglement generated between the 
self-trapping region and the propagative band by the ballis- 
tic motion of the quasi-free particles leaving the initial region. 
We also clearly see that in the SF-M case, the entanglement 
is vanishing in the Mott phase since there the state is a di- 
rect product of single occupancy local states. Notice also that 
the initial entanglement within the self-trapped condensate is 



basically conserved in time, showing a trivial time-evolution 
related to the periodic oscillations of the trapped condensate. 
Indeed, in figure [9] we have plotted the entanglement entropy 
profile at integer multiples of the Bloch period tb = 2tt / F 
and of the half-period to show the very nice superposition of 
the profiles for the trapped particles. Basically, it means that 
the initial correlations within the trapped wave-packet remain 
unchanged during the time evolution and no entanglement is 
lost or created. 
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FIG. 9. (Color online) Left: Entanglement entropy profile at different 
multiples of the Bloch period tb — 2ti/F and half-period tb/2 
for F = 0.3 with a SF initial state generated with A = 50 and 
/io = 0. Right: Entanglement entropy profile at different multiples 
of the Bloch period tb — 2ti / F and half-period rs/2 for F = 0.5 
with a SF-M initial state generated with A = 50 and /io = 2. 



IV. CONCLUSION 

We studied the dynamical behavior of a hard core boson gas 
initially prepared in a parabolic trap and subject to the sudden 
quench of a linear ramp potential just after the release of the 
trapping potential (like in the so-called Galileo ramp experi- 
ment with classical particles). The study was based on exact 
numerical diagonalization methods and on a hydrodynamical 
description which allowed the very good understanding of the 
dynamical behavior at the level of local density and local cur- 
rent density. The dynamics of this setup is rich, showing an 
escape of particles and Bloch oscillations for the self-trapped 
remaining condensate. Depending on the initial trap, the state 
of the condensate is either in a purely superfluid phase or in a 
mixed state with a Mott phase surrounded by two superfluid 
phases. In these two cases the behavior of the self-trapped 
condensate is different, showing in particular in the SF-M 
phase the possibility for a breathing condensate at high forces 
while in the SF situation the condensate is always oscillating 
as a whole. Notice finally that this setup can be used to cre- 
ate entangled many-body wave-packets by coupling the Bloch 
bound states to propagative modes. 

This work is supported by ANR-09-BLAN-0098-01. M. C. 
and P. W. benefited from the support of the International Grad- 
uate College on Statistical Physics and Complex Systems be- 
tween the universities of Lorraine and Leipzig. 



[1] F Bloch, Z. Phys 52, 555 (1928); C. Zener, Proc. R. Soc. Lon- [2] C. Waschke etal, Phys. Rev. Lett. 70, 3319 (1993); M. Ben Da- 
don Ser. A 145, 523 (1934). 



8 



han et al, Phys. Rev. Lett. 76, 4508 (1996); V. Agarwal et al, 
Phys. Rev. Lett. 92, 097401 (2004); O. Morsch a/. , Phys. Rev. 
Lett. 87, 140402 (2001); Z. He et aL, Phys. Rev. E 76, 056605 

(2007) . 

[3] L. D. Landau, Phys. Z. Sov^jetunion 2, 46 (1932); C. Zener, 

Proc. R. Soc. London 137, 696 (1932). 
[4] A. R. Kolovsky, Phys. Rev. Lett. 90, 213002 (2003). 
[5] M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, G. Rojas- 

Kopeinig, and H. C. Nagerl, Phys. Rev. Lett. 100, 080404 

(2008) . 

[6] L Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008). 

[7] J. Dziarmaga, Advances in Physics 59, 1063 (2010). 

[8] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, 

Rev. Mod. Phys. 83, 863 (2011). 
[9] F. Heidrich-Meisner, S. R. Manmana, M. Rigol, A. Muramatsu, 

A. E. Feiguin, and E. Dagotto,, Phys. Rev. A 80, 041603 (2009). 
[10] A. Minguzzi and D. M. Gangardt, Phys. Rev. Lett. 94 240404 

(2005). 



[1 1] X. Cai, S. Chen, and Y. Wang, Phys. Rev. A 84, 033605 (201 1). 
[12] M. Collura, H. Aufderheide, G. Roux, and D. Karevski, Phys. 

Rev. A 86, 013615 (2012). 
[13] D. Karevski, Eur. Phys. J. B 27, 147 (2002); G. M. Schuetz and 

S. Trimper, Europhys. Lett. 47, 164 (1999); S. Dorosz, T. Pla- 

tini, and D. Karevski, Phys. Rev. E 77, 05 1 120 (2008); T. Platini 

and D. Karevski, J. Phys. A: Math. Theor. 40, 1711 (2007); D. 

Karevski and T. Platini, Phys. Rev. Lett. 102, 207207 (2009); T. 

Platini, R. J. Harris, and D. Karevski, J. Phys. A: Math. Theor. 

43, 135003 (2010); Y. Ogata, Phys. Rev. E 66, 016135 (2002); 

66, 066123 (2002); W. H. Aschbacher and C.-A. PiUet, J. Stat. 

Phys. 112, 1153 (2003). 
[14] V. Hunyadi, Z. Racz and L. Sasvari, Phys. Rev. E 69, 066103 

(2004). 

[15] T. Platini and D. Karevski, Eur. Phys. J. B 48, 225 (2005). 
[16] R Calabrese and J. Cardy, J. Phys. A: Math. Theor. 42 504005 
(2009). 



